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hadronic Regge trajectories are essentially nonlinear and can be well approximated, 
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1 Introduction 



It is well known that the hadrons composed of light (n = (u,d),s) quarks populate 
approximately linear Regge trajectories; i.e., the orbital momentum £ of the state is 
proportional to its mass: £ = a'M 2 (£) + a(0), where the slope a' depends weakly on the 
flavor content of the states lying on the corresponding trajectory. Therefore, knowledge of 
Regge slopes and intercepts is very useful for spectral purposes. Since knowledge of Regge 
trajectories in the scattering region (t < 0) is also useful for many nonspectral purposes, 
e.g., in the recombination and fragmentation models, Regge trajectories become a valuable 
description of hadron dynamics, perhaps generally more important than the mass of any 
particular state. 

In the Veneziano model for scattering amplitudes [||] there are infinitely many excita- 
tions populating linear Regge trajectories. The same picture of linear trajectories arises 
from a linear confining potential P| and the string model of hadrons E9]. 
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However, the realistic Regge trajectories extracted from data are nonlinear. Indeed, 
the straight line which crosses the p and p^ squared masses corresponds to an intercept 
«p(0) = 0.48, whereas the physical intercept is located at 0.55, as discussed below. The 
nucleon Regge trajectory as extracted from the n + p backward scattering data is [|J 

a N (t) = -0.4 + 0.9£+- 0.25 1 2 , (1) 

and contains positive curvature. Recent UA8 analysis of the inclusive differential cross 
sections for the single-diffractive reactions pp — > pX, pp — > Xp at yfs = 630 GeV reveals 
a similar curvature of the Pomeron trajectory ||: 

a P (t) = 1.10 + 0.25* + - (0.16 ± 0.02)t 2 . (2) 

An essentially nonlinear a 2 trajectory was extracted in ref. || for the process n~p — > rjn. 
In addition to being disfavored by these experiments, linear trajectories also lead 
to problems in theory. Linear trajectories violate Cerulus-Martin fixed angle scattering 
bound 0]. (The so called square-root trajectory, used extensively in this paper for spec- 
troscopic purposes, saturates this bound ||.) Further, they violate the Froissart bound, 
in the following sense: The consequences of the strong duality between the saturation 
of the S'-matrix by i) narrow resonances || and ii) Regge asymptotic behavior JTIJ in 
the S'-matrix formulation of statistical mechanics by Dashen, Ma and Bernstein JTI[ was 
studied in [|I2" |. There it was shown that in order that the total cross section satisfy the 



Froissart bound |U|, a toi (s) < A\n s, A = const, the density of resonances in the elastic 



amplitude should decrease with energy, typically as ~ 1/E JL2|. (A decreasing density of 
resonances is a typical feature of dual amplitudes with Mandelstam analyticity (DAMA) 
||14|| .) A linear trajectory violates the requirement of a decreasing density of resonances 
(for a linear trajectory, the density of resonances grows as ~ E 3 [|HJ), and therefore, 
violates the Froissart bound in this framework. 

The idea of nonlinear Regge trajectories is not new. In the late-60's, by introduc- 
ing Regge cuts, through the eikonal method, a number of authors have shown that the 
effective trajectory for large momentum transfer (— t) goes like \f—i JTB|]. Subsequent 



comprehensive analysis by Vasavada, combining both Regge poles and cuts, arrived at 



the effective square-root trajectory in the entire complex angular-momentum plane |17 



The a(t) ~ y—i trajectory was also found by Gribov et al. |T8[ in a Pomeron- reggeon 



field theory. Other nonlinear forms have been studied also in the literature. 

Once the nonlinearity of Regge trajectories is an established fact, the determination 
of its actual form becomes an important issue. In this paper we address this issue. We 
start in Section 2 with a model study of a heavy quark- ant iquark system in a potential 
fitted to the results of an unquenched lattice QCD |19| , and show that its (parent) Regge 
trajectory is equally well fitted by both limiting cases of nonlinear forms allowed by 
dual amplitudes with Mandelstam analyticity (DAMA) |2(J (and thus, one can argue 
that all these nonlinear forms would fit the trajectory as well. The two limiting cases 
are the so called "square-root" trajectory and the "logarithmic" trajectory.) To gain 
some understanding of this peculiar observation, we turn to analytic model of a massless 



string with a variable tension (Section 3). (We show that this model is a relativistic 
generalization of a nonrelativistic rod with arbitrary potential, in the same sense as the 
usual Nambu-Goto string model is a generalization of a nonrelativistic linear potential 
model.) In this model we are able to recover underlying potential (or, in other words, the 
intergal of the variable string tension) from a known Regge trajectory, and we present the 
two potentials corresponding to the square-root and logarithmic trajectories, respectively, 
and compared them with the unquenched lattice QCD potential plf used in the heavy 



quark model. The two models (Section 2 vs. Section 3) are very different, but since the 



nonlinearity of Regge trajectories arises due to the flux tube breaking |19| , one can expect 
the same qualitative behaviour in both heavy and light quark systems. 

Sections 2 and 3 show that even though the two extreme nonlinear forms can both fit 
the bound states spectra, one of them, the sqaure-root form, goes beyond the position of 
the few lowest lying poles. This is why we choose the square-root form to fit and predict 
real-world spectra in Section 4. It should be noted that any nonlinear form bracketed by 
the square-root and logarithmic ones can be expected to give comparably good results 
for the lowest lying states, and a true test would be higher excited states (likely yet not 
observed) . 

With a few additional assumptions that are commented on as they are introduced, 
we obtain an excellent agreement with data. This means that the Regge trajectories are 
indeed nonlinear, and well approximated by the square-root form. We summarize and 
comment on some of the results at the end of Section 4. The last section contains our 
overall summary and conclusions. 

2 Heavy quarkonia model 



In a previous paper |T9| we considered a Hamiltonian model for heavy quarks which 
attempted to include the effect of light pair creation, or, in other words, the breaking of 
the flux tube stretched between the two heavy sources that inevitably occurs when the 
distance between the sources is sufficiently large. We have shown that consequently, Regge 
trajectories for the heavy quarkonia become nonlinear and their real parts terminate. In 
this paper we use the model to gain insight about the specifics of the (nonlinear) form of 
hadronic Regge trajectories. 

A meson consisting of two heavy quarks is well described in leading order by a non- 
relativistic, spin-independent Hamiltonian, viz. 

H = -i v2 + K(r) < (3) 

where m = M/2 is the reduced mass, with M = 5.2 GeV, and V(r) denotes a potential. 
We use M = 5.2 GeV, and a potential fit to results of unquenched lattice QCD calcula- 
tion for infinitely heavy sources. We have argued that our simple two body calculation 
produces results qualitatively similar to coupled channel method, because it dynamically 
takes into account coupling to open channels, and that the two approaches can be expected 
to differ significantly only for a few states near the threshold Jl9| . 



The screened static potential fitted to results of lattice calculations is [[ZT|]: 

V(r) = ( + ar) , 4 

V r J fir 

where [T x = (0.9±0.2) fm = (4.56±1.01) GeV" 1 , y/d = 400 MeV and a = 0.21±0.01. We 
obtain bottomonium-like spectra by diagonalizing the Hamiltonian (3) with the potential 
(4). The number of bound states in the model is finite, and it decreases with increased 
screening /i. The number of bound states is the largest for the parent trajectory, and 
decreases by one unit for each consecutive daughter. Another important feature of the 
Regge trajectories in this model is that they acquire curvature due to screening of the 
linear potential. Since we are interested in how the screening affects the form of the 
Regge trajectories (which would be linear in the absence of screening), for maximum 
reliability of the numerical fits it is useful to consider the case of a Regge trajectory 
with the largest number of bound states. Therefore, in what follows we concentrate on 
the parent trajectory for fi' 1 = 0.9 fm = 0.18 GeV -1 . The same conclusions as those 
presented below hold for all fi in the range indicated by the fit to the lattice data [P]J . 

A meson trajectory atjj(t), can be parametrized on the whole physical sheet in the 
following form: 

a fl (t) = a fl (0) + l\T'fi-(T fl -ty}, < v < \. (5) 



2 



(up to a power of logarithm), assuming that cxfi{t) is analytic function having a physical 
cut from to to oo, it is polynomially bounded on the whole physical sheet, and there exists 
a finite limit of the trajectory phase as \t\ — > oo |22j| . HereQ 7 is the universal slope which 
nonlinear trajectories must have in asymptopia |23|j , 

a(t) ~ —y(—t) u , \t\ — > 00; 

both 7 and the exponent v are flavor independent. cv,|(0) < 1, in accord with the Froissart 
bound. Note that for \t\ <C T, Eq. (5) reduces to the (quasi)linear form 

a fl {t) = a fl (0) + uiTYfH = a fl (0) + a' fi (0) t. (6) 

The subscripts i,j indicate which of the parameters depend on the quark content of 
the meson. In this section we drop the subscripts, since our model is applied here to a 
bottomonium-like system only. We consider the value of v restricted to lie between and 
1/2 (p^|, for details see Appendix A), v — should be understood as a limit v — > 0, ■yu 
fixed. In this limit, the difference of fractional powers reduces to a logarithm, viz., 

a (t) = a(0) - M log ( 1 - - J , 71/ = const. (7) 

Unlike a trajectory with v 7^ 0, the real part of the "logarithmic" trajectory does not 
freeze-out when t reaches T. The real part continues to grow, the only change for t > T 
is that the trajectory acquires a constant imaginary part. 

-'^Note that the choice of the signs in Eq. (5) is fixed by the requirements that a(t) must be real as 
t — > — 00, and have positive slope in the small-i region. 
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Figure 1: The parent trajectory for the bottomonium in our model compared to (a) our 
best fit of the "square- root" form, (b) our best fit of the "logarithmic" form. The solid 
vertical line corresponds to the threshold of the quenched lattice potential, i.e. (2M + 
a//i) 2 , the solid horisontal line shows the maximum £ for the data set. 
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The upper bound on v gives the so-called "square-root" trajectory, viz. 

a (t) = «(0) + 7[v / T- VT-t\. (8) 

When t reaches T, the real part of the "square-root" trajectory stops growing, and there 
are no states with a higher angular momentum than £ ma x = ot(T) . The parameter T is 
therefore the trajectory termination point. 

Fig. 1(a) shows the parent Regge trajectory for our model together with our best fit 
of the "square-root" form: 

a^(0) = -21.12 ±7.74, 

7^ = 2.68 ± 0.97 [GeV _1 ], 

T/- = 127.66 ± 1.98 [GeV 2 ] 

X 2 = 5.86 • 1(T 2 . (9) 

Our best fit of the logarithmic form, 



ai O g(0) 




= 


-10.89 ±7.59, 


(Ho, 




= 


5.0 ±3.92, 


T\ og 




= 


130.74 ±5.39 [GeV 2 




x 2 - 


= 1.45 -10~ 2 . 





(10) 

is shown in Fig. 1(b). In both cases, the large errors on the extracted parameters reflect 
the fact that there are only 8 data points for the trajectory. 

The two fits are indistinguishable in the region of the fit. In our case, since we know 
that in our model there is only a finite number of bound states, the square-root form 
(which terminates in general, and in our fit its termination point is near the termination 
of our bound-state data) is favored over the logarithmic form (which does not terminate). 
However, in practise, when one has no means to recognize whether the set of data is a 
part of a finite or an infinite ensemble, it is impossible to distinguish between the two 
forms based on fits to the bound state data. 

In view of this, we extend our considerations on how to numerically determine the form 
of a Regge trajectory from data. There is more information about Regge trajectories than 
just the masses of the bound states. Total cross sections can be related to the value of 
the relevant Regge trajectory at the origin (i.e., the intercept of the trajectory). In the 
region of negative t, the value of a(t) can be determined from the relevant differential 
cross sections. If there are reliable data at largeQ negative t in addition to bound states, 
it might be possible to obtain a more stringent constraint on the form of trajectory by 
fitting both sets of data simulataneously, providing the form is universally valid for all t. 



2 Unfortunately, according to our numerical simulations, data at larger negative t values, and with 
greater measurement precision, are both needed. 



3 Analytic model 

Our model calculation for heavy quarks revealed that the Regge trajectory formed by the 
bound states is certainly nonlinear but can be equally well fitted by a square-root form 
or a logarithmic form. 

In this section we try to shed some light on this observation using analytic model 
for a massless string with variable tension. The ends of the string can be massive. The 
model is intended to mimic a flux tube stretched between two quarks, and the varrying 
tension is to reflect possible dynamical effects such as weekening of the flux tube due to 
pair creation. 

The particular merit of the model for our purposes is that, given the form of a Regge 
trajectory, we are able to recover the form of underlying potential, in many cases in 
analytic form. We will show in this model why the two forms of Regge trajectories (i.e. 
square-root and logarithmic) are both likely to well approximate the Regge trajectory 
formed by the bound states of the unquenched lattice QCD potential (4). 

This section is organized as follows: first we review the standard relativistic Nambu- 
Goto string with massive ends and show that in the nonrelativistic limit this model 
simply corresponds to a rigid rod with a linear potential generated between its massive 
ends. After introducing the generalized string model, we show that it corresponds in the 
nonrelativistic limit to a rigid rod with arbitrary potential. Details fo the dynamics of 
the generalized string with massive ends, such as derivation of expressions for its energy 
and orbital angular momentum, can be found in the Appendix B. The key part of this 
section is devoted to the generalized string with massless ends. Within this framework we 
find potentials that lead to square-root and logarithmic trajectories, respectively. Since 
the dynamics of mesons (such as the termination of Regge trajectories) is dominated by 
the behavior of the flux tube, we claim that what we learn from the string with massless 
ends is qualitatively relevant to the case of massive quarks as well. 

3.1 The Nambu-Goto string and the generalized string. 

The action of the standard relativistic Nambu-Goto string with massive ends in the 
parametrization r = t = x° is written as (a is the string tension) [[| 

S = -a f 2 dt r ds^Jx' 2 (l - x 2 ) + {xx'f - J2 m i [ ' dt^fl- i?, (11) 

Jti JO £_i 2 Jtl 

X = X = x(t,s), Xi=x(t,Si), i = 1,2, Si = 0, S 2 = 7T, 

where from now on the dot and the prime stand for the derivative with respect to t and 
s, respectively, unless otherwise specified. 

To illuminate the physical meaning of this model, let us consider for the moment its 
nonrelativistic limit. In the nonrelativistic limit, \x(t, s)\ <C 1, |xj| <C 1, and the action 
reduces to || 

S = -a dt dsVx^- VmJ dt + V — / dtx 2 . (12) 

Jt! JO i=12 Jt X i=l2 2 J tl 



Integration over s gives the length of the string (under the assumption that there are no 
singularities on the string). The variation of the first term in the nonrelativistic action, 
Eq. (11), with respect to the string coordinates leads to the requirement on the string to 
have the form of a linear rod connecting the massive ends. The effective action that leads 
to the equations of motion of the massive ends is therefore 

Seff = £dt (-cr\ Xl (t) -x 2 (t)\ + E ^ ] • (13) 

Hence, in the nonrelativistic limit, the string generates a linear potential between its 
massive ends: V(\x\ — X2I) = cr\%\ — %t\- 

Here we generalize the standard string formulation described above to the case of 
an arbitrary potential between the string massive ends. Such generalization is done by 
the modification of the standard, constant, string tension into an effective string tension 
which is a function of |x|, as follows: 

S gen = - J* dt r dsa{\x\)Jx' 2 (l - x 2 ) + (xx') 2 - V m % [ * dtJl-x 2 . (14) 

Jti Jo i=l2 Jti 

The action is similar to that of the standard Nambu-Goto string. In the nonrelativistic 
limit, however, in place of (13) one will now obtain 



Vn,e// = j t dt [~ j dsa(\x\)\x'\ + E 



m,iX\ 

i=l,2 



t2 dW-y ( | Xl -x 2 |)+£^i]. (15) 



U * t=l,2 



In contrast to the previous case of the standard string, it is seen in the above relations 
that now 

U U \x'\ ds d\x\ ' l j 

i.e., the effective string tension is the derivative of a potential with respect to the distance. 
Obviously, in the case of a linear potential, the effective string tension reduces to the 
standard (constant) one. 

Similarly to the standard case of a constant string tension which represents the rela- 
tivization a la Poincare of the nonrelativistic two-body problem with linear potential f25fl , 
the generalized string can be considered as the relativization of a nonrelativistic two-body 
problem with an arbitrary potential. Details of the dynamics of the generalized string 
model are given in Appendix B. 



3.2 Generalized massless string 

Since we are interested in dynamical issues related to the properties of the flux tube, 
without a loss of generality it is sufficient for us to consider the case of a generalized 
string with a massless ends. This case is obviously more tractable than the case with 
arbitrary massess, and nevertheless exhibits the same qualitative features. 

The energy and orbital momentum of the generalized massless string, mi = rri2 = 0, 
is given by 



R d P a (p] J = 2 [ R d P a ^ UJ P 2 

y/1 — UJ 2 p 2 ' JO y/1 — LU 2 p' 2 



E = *l. -P^h. J = H ^P=&. < 17 > 



where R = 1/u) is half of the string length for a given u. The condition uR = 1 follows 
from, e.g., Eqs. (B.14) with rrii — > 0. 

By eliminating uj from Eqs. (17) one can obtain J as a function of E 2 , the Regge 



trajectory. It will be shown elsewhere [261 that it is possible to uniquely recover the 



potential (V(r) ~ / dp <j(p)) from the known analytic form of Regge trajectory, for both 
massless and massive generalized strings, and the techniques of the corresponding inverse 
problem will be presented in detail. In this way, we have recovered the potentials that 
correspond to both the square-root and logarithmic Regge trajectories. For our present 
purposes, here we present only the final results. 

3.2.1 Square-root trajectory 

The potential (a, /i = const, V(p) — > a/2p as p — > oo, and hence E — > a / p) 

V(p) = — arctan(7r/ip), (18) 

Tip 



for which 



dV(p) _ o 

dp 1 + ijipp)' 



"W = -^ = ^^W. (W) 



leads, via (17), to 



e= , Z , „ J = —Ji- . " „ ■ (20) 

y'U 2 + 7l 2 p 2 Tip 2 \ \JUJ 2 + TT 2 p 2 I 



Eliminating uj from the above relations gives 

1 



J = — [a/p-^/(a/p) 2 -E 2 ), (21) 

i.e., the square- root Regge trajectory. For E <ti o~/p, it reduces to an (approximate) linear 
trajectory, J ~ E 2 /(27ra). 
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Figure 2: Potentials which in our analytical model lead to logarithmic and square-root 
Regge trajectories, respectively, compared to the unquenched lattice QCD potential used 
in our heavy-quarkonia calculation. 

3.2.2 Logarithmic trajectory 

The potential 



tw n cr / , . logfl + (27iap) 2 } 

V{p) = 2 arctan(27r/ip) - - L 



2lTfI 



2npp 



(22) 



for which 



a(p) = a 



log[l + (2Tippf 

(27TPP) 2 



leads to 



E 



a 



. ., ,, > T a uj + \Ju 2 + 47r 2 /i 2 
u 2 + 47f 2 p 2 -u\, J=^—, 2 lo & 



2np 2 V J 27r/i 2 ° 2uj 

from which eliminating u (viz., uo = tt(<j 2 — p 2 E 2 )/(aE)) gives 

E 2 \ 



(23) 



(24) 



J 



2vr/i 2 



log 1 



(<r/»< 



(25) 



i.e., the logarithmic Regge trajectory. For E <C cr/p, it again reduces to an (approximate) 
linear form, J ~ E 2 j{2-no\ 
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Fig. 2 shows both potentials together with the confining part of the screened lattice 
QCD potential that we used in the heavy quark model in the previous section. All 
potentials are normalised to the same asymptotic value. The vertical lines show the 
region of the distance r which is relevant for the bound states. 

All three potentials are very close in the region of bound state physics. Therefore, 
one can conclude that they would lead to similar Regge trajectories as far as the bound 
states not close to the trajectory thresholds are concerned. One could further speculate 
that since the lattice potential is steeper than the square-root potential (which in turn, is 
steeper than the logarithmic potential) in asymptotia, it is reasonable to expect that the 
bound states of the lattice potential terminate at a lower t than the square-root trajectory, 
and that the square-root trajectory approximates the lattice potential trajectory better 
for states near the trajectory threshold than the logarithmic trajectory. (Recall that the 
real part of the logarithmic trajectory does not terminate.) 

Combining the conclusions of the analytic model for the generalized string with those 
of the heavy quark model of the previous section, we conclude that even though both non- 
linear forms considered may appear to be a good approximation to QCD, the square-root 
form is likely to be more realistic. Therefore, we use the square-root form for phenomeno- 
logical purposes. 

4 Trajectory parameters and spectroscopy 

In this section we determine, assuming that Regge trajectories are of the form (8), tra- 
jectories thresholds and intercepts using various experimental information. Typically, we 
use masses of a few known lowest lying states, and in the case the p trajectory we also 
use the value of the intercept (which is known and well-established) found from exchange 
processes. The value of 7 (the universal asymptotic slope) is fit to the p trajectory, and 
then taken as universal for all other trajectories. 

Whenever possible, we try to use as inputs states that are believed to be pure quark- 
antiquark states with definite flavors. For example, we do not use the masses of the 
observed 0, / or 77(77') states to find the parameters of ss trajectory because the physical 
states are mixtures of ss with the light quark-antiquark components. 

In some cases, e.g. for tensor mesons and axial-vector mesons, there are not enough 
data to determine the parameters. We are forced to make aditional assumptions. In 
particular, we assume that the thresholds of parity partner trajectories coincide. We 
comment further on this assumption where it is introduced in Section 4.2. Here we wish 
to mention that our additional assumptions are justified a posteriori by our excellent 
results. 

The approach has more predictive power than one would naively expect. This is 
because the parameters for different flavors are related by 

(i) additivity of intercepts, 

a {i (0) + a j - j (0) = 2a jl (0), (26) 
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where i,j(i,j) refer to the quark (antiquark) flavor^, and 
(ii) additivity of inverse slopes near the origin, 

-^ = 4-, (27) 

a- a- a- 

which is favored over another constraint suggested in the literature, factorization of slopes, 



J3 v Ji> 



a il a 'il = ( a 'il) i ^°y ^ ne neav y quark limit [p7|1 . These two additivity requirements are 



independent of which specific form is assumed for the trajectories. 

Once the parameters of the square-root trajectories are known, the first obvious appli- 
cation is meson spectroscopy. We calculate masses of a few excited states lying on each of 
the trajectories under consideration, and compare with experimental data from ref. |28j . 
Where there are no data available, our results are predictions. 

We start with vector meson trajectories, and then consider tensor meson, pseudoscalar 
meson, and finally, axial-vector meson trajectories. The section concludes with a brief 
discussion of the results. 

4.1 Vector mesons 

We start with the p trajectory. The intercept of this trajectory is well established. From 
the behavior of the differential cross section of the process n~p — > n n the intercept a p (0) 
has been extracted to be 0.58 ||, 0.56 ± 0.01 0, 0.56 ± 0.02 g|, 0.53 [||. Extensive 



analysis of this process by Hohler et al. |33| leads to a p (0) = 0.55 [B3]. From the 



difference of the total cross sections of 7t + p and tt p scattering a p (0) has been inferred to 
be 0.57 ±0.01 [55| , 0.55 ±0.03 |36]]. Bouquet finds from the dual topological unitarization 



[37] 0.51 < Q!p(0) < 0.54. All of the above values are consistent with 

a p (0) = 0.55, (28) 

in agreement with the value of the p trajectory intercept extracted by Donnachie and 

Landshoff from the analysis of pp and pp scattering data in a simple pole exchange model 

3~8| . We therefore take the value of the intercept, Eq. (28) as one of the constraints on 



the form of the p trajectory. Two more constraints that are needed are provided by the 
mass and spin of p and p$, i.e.: (i) a p {M p v ) = 1, (ii) a p (M p 2 a ) = 3, where M p = 769.0 ±0.9 



MeV and M P3 = 1688.8 ±2.1 MeV 



Inserting these values into the functional form (8), we extract the parameter values 
7 = 3.65±0.05GeV-\ J¥ p = 2.46 ± 0.03 GeV. (29) 



It is interesting to note that the values of parameters we find based on spectroscopy and 
the known intercept are in excellent agreement with 



7 



3.72 ±0.30 GeV" 1 , Jt p = 2.50 ± 0.10 GeV, (30) 



3 This additivity is a firmly established theoretical constraint on Regge trajectories (for an extensive 
list of references see [p7[). 
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extracted in ref. JJjJ for a similar form of the p trajectory from the analysis of ttN 

charge-exchange scattering data.Q Earlier analysis of ref. [ID| found JT p = 2.4 ±0.4 GeV. 

The parameters of the K* trajectory are obtained by using ax*{Mj c ,) = 1 with 

M K *o = 896.1 ± 0.3 MeV and a K ,(M 2 K .) = 3 with M K *o = 1776 ± 7 MeV p8) , and 



taking the value of 7 found from p spectra as a universal slope in asymptopia. This yields 



T K * = 2.58 ± 0.03 GeV. 
a K *(0) =0.414 ±0.006 



(31) 



The value of the intercept is in an excellent agreement with the results of the analysis of 
hypercharge exchange processes 7i + p — > K + T 1 + and K~p — > 7r~S + | |41| |. 

Taking the parameters of the p and K* trajectories as known, those for the trajectory 
may be obtained from the requirements of additivity of inverse slopes and intercepts, Eqs. 
(26),(27): 



tty(0) =0.25 ±0.02, 
fTH, = 2.59 ± 0.11 GeV. 



(32) 



Similarly, one can obtain the parameters of trajectories for the states containing c- 
and 6-quarks, using the value of the universal slope 7 found from p-trajectory and masses 
of the corresponding states. In particular, we use masses of D*, J ftp, B* and T as inputs. 
The remaining trajectories are determined from the requirements of additivity of inverse 
slopes and of additivity of intercepts, Eqs. (26), (27). The parameters of these vector 
meson trajectories are summarized in Table I. 





P 


K* 







a(0) 


0.55 


0.414 ±0.006 


0.27 ±0.01 




y/f, GeV 


2.46 ±0.03 


2.58 ±0.03 


2.70 ±0.07 






D* 


D* s 


J/1> 




a(0) 


-1.02 ±0.05 


-1.16 ±0.05 


-2.60 ±0.10 




y/T, GeV 


3.91 ±0.02 


4.03 ±0.04 


5.36 ±0.05 






B* 


B* s 


B* 


T 


a(0) 


-7.13 ±0.17 


-7.27 ±0.17 


-8.70 ±0.18 


-14.81 ±0.35 


\/T, GeV 


7.48 ±0.02 


7.60 ±0.04 


8.93 ±0.03 


12.50 ±0.02 



Table I. Parameters of the vector meson trajectories of the form (8). The intercept of 
the p trajectory was taken as an input. 

Using the parameters shown in Table I, we calculate masses of the spin-1, spin-3 and 
spin-5 states lying on these trajectories. Our results are compared with data from [^] 
in Table II. In this, and subsequent Tables IV, VI, VIII, the values used as input for our 
analysis are shown in boldface. The masses of the states for which there are no data 
available should be considered as our predictions. 



4 The p trajectory adopted in J39|] contains an additional term — O.li^iAI^ — t 1 to take into account 
nonzero resonance widths, which reduces the planar intercept a p (0) down to 0.51. 
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J= 1 


J = 


= 3 


J = 5 




This work 


ref. |8| 


This work 


ref. |8| 


This work ref. [28] 


Oi p (t) 


769.0 ± 0.9 


769.0 ±0.9 


1688.8 ±2.1 


1688.8 ±2.1 


2124 ±19 


a K *{t) 


896.1 ±0.3 


896.1 ±0.3 


1776 ± 7 


1776 ± 7 


2215 ±21 


a<t,{t) 


1015 ±17 


1019.4 


1863 ±31 


1854 ± 7 


2305 ± 42 


a D *{t) 


2006.7 ±0.5 


2006.7 ±0.5 


2721 ± 23 


3191 ±22 


a D * (t) 


2102 ±29 


2106.6 ±2.1 ±2.7 


2808 ± 28 


3279 ± 30 


<Xj/i/>(t) 


3096.9 


3096.9 


3753 ± 41 


4240 ± 39 


a B *(t) 


5324.9 ±1.8 


5324.9 ±1.8 


5814 ±51 


6217 ±46 


a B ; (t) 


5411 ±58 


5416.3 ±3.3 


5901 ±53 


6306 ± 49 


a B * (t) 


6356 ± 80 


6853 ± 72 


7276 ± 65 


a?(t) 


9460.4 ± 0.2 


9460.4 ±0.2 


9906 ±91 


10304 ±84 



Table II. Comparison of the masses of the spin-1, spin-3 and spin-5 states given by ten 
vector meson trajectories of the form (8) with data. All masses are in MeV. 



4.2 Tensor meson trajectories 

Tensor meson trajectories are parity partners of the vector meson trajectories. In a nonrel- 
ativistic theory, tensor meson trajectories are degenerate with vector meson trajectories. 
In a field theory, both acquire different corrections to the intercept, but they can still be 
expected nearly degenerate in the bound state region. From the form of the trajectory, 
Eq. (8), and since the intercept is in practise much smaller than the other terms in Eq. 
(8), it is clear that trajectories can be near-degenerate in the bound-state region even if 
their intercepts differ by a large percentage, but not if their thresholds are very differ- 
ent. For this reason, and to reduce the number of free parameters, we assume that the 
thresholds of tensor meson trajectories are the same as the thresholds of the vector meson 
ones. 



TV 



Our assumption is supported by analysis of the processes n + p 



K+E+ and K~ 



P 



41 1 . Most of the analysed experiments show the same slope at the origin for K* 



and K\ trajectories. In our formalism, this translates to equal thresholds. 

To further test this assumption, we first fit the threshold of the K% trajectory which is 
the only one with more than one state. (A recent measurement of the a^ mass [33, [ 
four standard deviations below the previous value fl28|, see discussion in the text below. ) 



is 



The parameters of the K\ trajectory are fixed by using 7 determined for p trajectory, 
ajc*(M£.) = 2 with M K *o = 1432.4 ± 1.0 MeV and a K *(Ml*) = 2 with M K *o = 2045 ± 9 
MeV [EH]. The calculation yields 



T, 



A'.f 



2.64 ± 0.03 GeV 



(33) 



The value of the K% trajectory threshold is consistent with that of the K* trajectory given 
in Table I, supporting our assumption. Moreover, it is in excellent agreement with the 
values for the slope shown in ref. ||l] . 
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We therefore conclude that our assumption that the thresholds of the parity partners 
are the same is plausible. Under this simplifying assumption only one state on each 
trajectory is needed to completely fix the parameters. We choose to fit a 2 , K%, Xc2 and 
Xb2 trajectories using the masses of the corresponding lowest lying state. The remaining 
trajectories are determined using the addititivity requirements. Parameters found in this 
way are presented in Table III. 





a 2 


K* 2 


/ 2 




o(0) 


0.60 ±0.03 


0.42 ±0.03 


0.22 ±0.07 






D$ 


n* 


Xc2(lP) 




a(0) 


-1.16 ±0.05 


-1.35 ±0.05 


-2.93 ±0.09 






B* 2 


B* s2 


B* c2 


Xb2(lP) 


a(0) 


-7.62 ±0.13 


-7.81 ±0.13 


-9.39 ±0.14 


-15.85 ±0.26 



Table III. Parameters of tensor meson trajectories of the form (8). 
As before, we use the knowledge of Regge trajectories for spectroscopy purposes. Our 



results are compared with data from E3] in Table IV. The masses of the states for which 



there are no data available should again be considered our predictions. 





J = 


--2 


J = 


--A 


J- 


= 6 




This work 


ref. @ 


This work 


ref. H| 


This work 


ref. |2| 


a a2 (£) 


1318.1 ±0.6 


1318.1 ±0.6 


1927 ±18 


2020 ± 16 


2256 ± 21 


2450 ± 130 


a K ;(t) 


1432.4 ±1.3 


1432.4 ±1.3 


2026 ± 20 


2045 ± 9 


2357 ±24 


a ti(t) 


1544 ± 37 


1525 ±5 


2124 ±40 


2457 ± 48 


old* (t) 


2454 ± 23 


2458.9 ±2.0 


3010 ±22 


3390 ± 21 


a D *Jt) 


2560 ±27 


2573.5 ±1.7 


3109 ±29 


3489 ± 31 


a Xc2(lP)( t ) 


3556.2 ±0.1 


3556.2 ±0.1 


4092 ± 39 


4498 ± 38 


dB* (t) 


5698 ± 45 


5698 ± 12 


6122 ±42 


6472 ± 39 


<*B*Jt) 


5797 ±47 


6220 ± 45 


6570 ± 43 


(*B*Jt) 


6780 ± 52 


7213 ±50 


7582 ± 49 


° ! Xb2{lP)( t ) 


9913.2 ±0.6 


9913.2 ±0.6 


10310 ±72 


10665 ± 68 



Table IV. Comparison of the masses of the spin-2, spin-4 and spin-6 states given by ten 
tensor meson trajectories of the form (8) with data. All masses are in MeV. See the text 
regarding the a 4 and a 6 masses. 

Note that all the masses agree very well for the available data, except for M a4 . With 
respect to this, we wish to mention that recent analyses by the VES collaboration of the 
reactions ir~Be — > 7f + 27T~27r Be [42] and ir~Be — > 7i + 27i^2Be ^j reveal the mass of the 
a 4 resonance as seen in the a 4 — > up and a 4 —* f 2 7i decays, respectively: 1944 ± 8 ± 50 



MeV ||2f and 1950 ± 20 MeV 0. Our value in Table IV is in excellent agreement with 
these latest measurements. 

According to our fitted square-root form for the a 2 trajectory, the last state below 
threshold is J = 8. This conclusion is very sensitive to the functional form assumed, even 
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though the threshold value itself is not (see Section 2). Generally, as the mass approaches 
the threshold, we expect to find larger discrepancies between the observed states and 
predictions of any specific form of trajectory, Eq. (5). Nonetheless, the discrepancy 
for the next-to-last, J = 6 state appears to be less than 10%. Note also that our K\ 
prediction appears to be below the data, consistent with a^ and a 6 deviations, suggesting 
a possible systematic effect. This raises the possibility that the growth of the square-root 
form of trajectory near the termination point is too rapid and that the true trajectories 
are somewhat flatter in this region. Thus, it is possible that the J = 6 state is actually 
the last state on these trajectories. This point of view is further supported by the large 
width of the a&, as would be expected for a state near the termination point. Hence we 
also expect the K$ to be similarly broad. 

4.3 Pseudoscalar meson trajectories 

Here we calculate the parameters of the pseudoscalar meson trajectories. We start with 
the pion trajectory and fix its intercept and threshold by using the masses of the two lowest 
lying states, M w o = 135 MeV and M n2 = 1677 ± 8 MeV @, in relations a n (M%) = 0, 
a^{Ml 2 ) = 2, and 7 from Eq. (29). 

Similarly, the intercept and the threshold of the iT-trajectory are determined using 
the massses of K and K 2 . 

The parameters of charmed mesons trajectories are found using the masses of D and 
rj c , utilizing the additivity requirements, Eqs. (26), (27), and the parameters of the light 
mesons (7r-trajectory) found above. 

If one tries to procceed similarly in the case of mesons containing the b quark, and use 
the only two experimentally known masses as inputs (the mass of B and B s ) in conjunction 
with the parameters of n- and .^-trajectory, the error estimates on the extracted values 
exceed two hundred percent. The results turn out to be highly unstable with respect to 
the threshold of the 7r-trajectory. This is a peculiar consequence of the pion intercept 
being close to zero, and J = of the lowest lying state. 

Conversely, this observation allows us to convert the problem into an additional self- 
consistency check. We use the 6-quark meson spectra to extract information about the 
it trajectory To do that, an additional input is needed. We use the mass of rjb given 
in ref. ||44j| , and fit the value of the pion threshold in addition to the parameters of 
^-containing mesons. Recall that the light meson parameters affect the heavy quark 
parameters only through the additivity requirements. Since the additivity requirements 
are well established, such a self-consistency check is quite nontrivial. 

The parameters of the pseudoscalar trajectories that we obtain are shown in Table V. 
The parameters of the 7r-trajectory given in the Table come from fitting the pion spectra. 
From the bottom spectra we obtain the following value of the pion threshold: 



'T w = 2.79 ± 0.40 GeV (34) 

Note the agreement with the value given in Table V. 
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7T 


K 


Vs 




a(0) 


-0.0118 ±0.0001 


-0.151 ±0.001 


-0.291 ±0.003 




VT, GeV 


2.82 ±0.05 


2.96 ±0.05 


3.10 ±0.11 






D 


D s 


Vc 




a(0) 


-1.61105 ±0.00005 


-1.751 ±0.001 


-3.2103 ±0.0001 




VT, GeV 


4.16 ±0.03 


4.29 ±0.06 


5.49 ±0.02 






B 


B s 


B c 


Vc 


a(0) 


-7.41 ±0.17 


-7.54 ±0.17 


9.00 ±0.17 


-14.80 ±0.34 


VT, GeV 


7.89 ±0.16 


8.01 ±0.16 


9.24 ±0.12 


12.98 ±0.24 



Table V. Parameters of the pseudoscalar meson trajectories of the form (8).[| 



Our results for pseudoscalar spectroscopy are compared with data from |28[ in Table 
VI. The masses of the states for which there are no data available should once again be 
considered our predictions. 





J = 


J = 2 


J = 4 




This work ref. [2£] 


This work ref. [2~B| 


This work ref. |2~B 


a w (t) 


135 135 


1677 ±8 1677 ±8 


2237 ±26 


a K (t) 


493.7 493.7 


1773 ± 8 1773 ± 8 


2333 ± 27 


a ris 


698 ± 14 


1869 ± 38 1854 ± 20 


2429 ± 54 


a D (t) 


1864.1 ±1.0 1864.1 ±1.0 


2692 ± 19 


3228 ± 22 


a Ds (t) 


1971 ±19 1969.0 ±1.4 


2786 ± 26 


3323 ± 32 


(*vM 


2979.8 ±2.1 2979.8 ±2.1 


3692 ± 23 


4217 ±25 


a B {t) 


5279.8 ±1.6 5279.8 ±1.6 


5830 ± 89 


6286 ± 93 


®bM 


5369.6 ±2.4 5369.6 ± 2.4 


5920 ± 89 


6376 ± 93 


a Bc {t) 


6283 ± 79 


6826 ± 79 


7287 ±80 


a Vb (t) 


9424 ±3.6 


9914 ± 148 


10353 ±150 



Table VI. Comparison of the masses of the spin-0, spin-2 and spin-4 states given by ten 
pseudoscalar meson trajectories of the form (8) with data.0 All masses are in MeV. 



4.4 Axial-vector meson trajectories 

Finally, we calculate the parameters of the axial-vector trajectories. These have the same 
parity and spin as the parity partners of the pseudoscalar trajectories - the pseudovector 
trajectories^. The available data are insufficient to fix both thresholds and intercepts of 
axial-vector trajectories. We assume, therefore, in analogy with the tensor vs. vector 

5 The parameters for the -ftf-trajectory were found using the mass of K^ from [£8|. If we instead use a 
mass of the corresponding pure ns state as found in ref. [|45| , i.e. Mr- 2 — 1762 ± 18 GeV, the parameters 
change slightly: the intercept -0.153 ± 0.003, and the threshold 2.93 ± 0.07 GeV. 

6 We take the error estimate on the 775 mass as 10% of the calculated splitting, in agreement with Fig. 
2 of the second paper of ref. [Q . 

7 We use the usual Regge terminology for parity partners, which is not the same as normally used in 
lattice QCD/chiral condensates. 
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meson trajectories, that the thresholds of the parity partners coincide, and further, that 
the thresholds do not depend on charge conjugation in accordance with the C-invariance of 
QCD. Under this assumption, we use the calculated thresholds of pseudoscalar trajectories 
and fit only the intercepts of the axial- vector trajectories using the masses of a few lowest 
lying states and/or the additivity requirements, Eqs. (26), (27). 

In particular, we use masses of Oi, D s i, x cl (lP) and x fel (lP). Note, that to find 
the parameters of trajectories of mesons containing strangeness, we use in the absence 
of data for K\ the mass of D s rather than f v This is because the f 1 is not a pure ss 
state, and therefore the value of the intercept extracted from the physical state could 
deviate significantly from the value corresponding to a pure ss trajectory. (An additional 
complication arises, in principle, for axial mesons since the physical qq' states are mixtures 
of axial-vector and pseudovector components, because charge conjugation is not well- 
defined for other than qq states. Note that for all states other than axial-vector mesons, 
the corresponding state with opposite C is exotic, so that the problem is avoided. Of the 
inputs we use, this is relevant only to -D sl . We assume that the masses of D s ia and D s i B 
are close, because the heavy quark content implies that the hyperfine mass splitting is 
small.) 

As before, 7 is taken as universal, and the value found from the p trajectory is used. 
The parameters we obtain in this way are in Table VII. 





a\ 


Kt 


A 




a(0) 


-0.03 ±0.07 


-0.22 ±0.08 


-0.42 ±0.15 






D x 


D sl 


Xd(lP) 




a(0) 


-1.83 ±0.05 


-2.03 ±0.06 


-3.63 ±0.07 






B{ 


B sl 


B c \ 


XbiO-P) 


a(0) 


-7.87 ±0.27 


-8.06 ±0.28 


-9.67 ±0.27 


-15.70 ±0.54 



Table VII. Parameters of axial- vector meson trajectories of the form (8). 
Our results for masses of the states lying on the resulting trajectories are compared 



with data from 28 in Table VIII. The masses of the states for which there are no data 



available should, as before, be considered as our predictions. 





J = 


= 1 


J = 3 


J = 5 




This work 


ref. HH 


This work ref. |f28 


This work ref. 


m 


oc ai (t) 


1230 ± 40 


1230 ± 40 


2000 ± 31 


2427 ± 32 


a KlA (t) 


1368 ± 46 


2109 ±33 


2535 ± 33 


a h {t) 


1501 ± 78 


1518 ±5 


2218 ±62 


2643 ± 67 


a Dl (t) 


2418 ±26 


2422.2 ±1.8 


3042 ± 25 


3473 ± 24 


OCD sl (t) 


2535.4 ±0.3 


2535.4 ±0.3 


3150 ±34 


3580 ± 37 


a Xcl(lP)(*) 


3510.5 ±0.1 


3510.5 ±0.1 


4080 ± 30 


4513 ± 28 


a Bl (t) 


5692 ± 104 


6171 ±102 


6570 ± 104 


«B.i(*) 


5796 ± 105 


6273 ± 104 


6671 ± 105 


a Bcl (t) 


6740 ± 95 


7215 ± 92 


7620 ± 90 


°hcbi(iP)(t) 


9891.9 ±0.7 


9891.9 ±0.7 


10333 ± 173 


10727 ±172 
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Table VIII. Comparison of the masses of the spin-1, spin-3 and spin-5 states given by 
ten axial-vector meson trajectories of the form (8) with data. All masses are in MeV. 



4.5 Discussion of the results 

Here we comment on the results presented in this section: 

4.5.1 Intercepts 

Light- quark (nn) trajectories: 

i) The intercept of the linear trajectory that passes through p and ps, a(t) = 0.48 + 0.88£, 
should be reduced for the a 2 if the a 2 were on a parallel linear trajectory: 0.47 + 0.88£. 
Even lower value for the intercept (and a different value of the slope), 0.30 + 0.98 t, is 
needed if a 2 and a^ masses were to lie on the same linear trajectory (the latter being 1.944 
MeV, in accord with VES's recent measurements [0). (Note the discrepancy between the 
two sets of the a 2 trajectory parameters.) In our case, due to nonlinearity the intercept 
increases to 0.55 for p and 0.60 ± 0.03 for a 2 . The data on pp, pp [f5J undoubtedly show 
that the tensor intercept is larger than the vector one (although the numerical errors 
of these data prevent us from using them as input). Therefore, data disfavor a linear 
trajectory, 
ii) The value a a2 (0) = 0.60 ± 0.03 that we obtain here is in agreement with the recent 



reanalysis by Cudell, Kang and Kim [^6] of the pp and pp scattering data in the simple 
pole exchange model |38| « p (0) = 0.50 ± 0.07, a a2 (0) = 0.66 ± 0.08. In this analysis the 
pairs of trajectories, p and u>, and a 2 and f 2 , are taken to be degenerate. However, since 
M(uj) > M(p) and M(f 2 ) > M(a 2 ), the intercepts of these trajectories should actually 
satisfy a w (0) < a p (0) and a a2 (0) < a/ 2 (0). Hence, the value of a p (0) should be larger the 
average of a p (0) and a^(0), i.e., 0.50. Similarly, the value of a a2 (0) should be less than 
0.66. 

Our values for the intercepts agree well with their values within the errors, moreover, 
they satisfy the above described inequalities, 
iii) Our value a ai (0) = —0.03 ± 0.07 is consistent with a ai (0) = —0.1 ± 0.2 found in the 



analysis of the reaction n~p — > n^n+n |47| . 

Trajectories containing strangeness: 
The values of K* and K 2 intercepts are in excellent agreement with the analysis of hy- 
percharge exchange processes 7r + p — > K + H + and K~p — > 7r~£ + p4 . 



Trajectories containing a heavy quark: 
i) The calculated intercepts of the heavy quark trajectories are negative, and the absolute 
value of an intercept increases with the increasing mass of the heavy quark. This feature 
is confirmed by means of the operator expansion of QCD and dispersion relations for 
heavy quark four-current correlation functions by Oganesyan and Khodzhamiryan p8| . 
They also find the estimates aj/^(0) ~ —(2 — 3), q;t(0) < —10. Our corresponding results, 
a J/^(0) = —2.60 ± 0.1 and «t(0) = —14.81 ± 0.35, agree very well with these estimates 
provided by QCD. 
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ii) Our predictions for the intercepts of the charmed meson trajectories may be confronted 
with the existing data. For the reaction 

it- A -»• D(D*)X 

the differential cross section, as given in the triple- Regge limit [IS], is 

d 2 a 



dxFdp\ 



A{l-\x F \) n e- bp ±, A, b = const, 



where xp and p± are the Feynman-x and transverse momentum variables, respectively, 
and n = 1 — 2a(0), with a(0) being the intercept of the trajectory exchanged in the 
t-channel of the reaction. The D meson production proceeds via the exchange of both D* 
and D\ trajectories. Hence, on the basis of our results for the corresponding intercepts, 
one may expect n ~ 3.2 ± 0.2 for the -D-production. In 7r~H interactions at 360 GeV 



the NA27 experiment by the LEBC-EHS collaboration Q finds n = 3.80 ± 0.63. For 



350 GeV ix beam on emulsion the WA75 collaboration [51] finds n = 3.5 ± 0.5. In tx Cu 



interactions at 230 GeV the ACCMOR collaboration J52| finds n = 3.23 ±0.29 for leading 
charm production (and n = 4.34 ± 0.35 for nonleading, and n = 3.74 ± 0.23 for the 
combined data). For 340 GeV ix~ beam on Si, Cu and W targets the WA82 collaboration 
fx| finds n = 2.9 ± 0.3. Finally, for 250 GeV ir~ beam on Be, Al, Cu and W targets the 
E769 collaboration [Q finds n = 3.2 ± 0.5 on Al target, and n = 3.7 ± 0.4 for leading 
D production on all targets (and n = 4.0 ± 0.4 for nonleading, and n = 3.9 ± 0.3 for the 
combined data). 

Since the D* meson production can proceed via the exchanges of all four, D*, D^, D 
and Di, trajectories calculated here, definite predictions for n are more complicated than 
in the previous case. Of course, quantitatively n should be larger in this case, because 
the intercepts of D and Di trajectories are larger in magnitude than those of D* and D\ 
trajectories. In fact, in the experiment mentioned above the E769 collaboration ]S3| finds 
n = 3.5 ± 0.3. Earlier measurements in the NA27 experiment []50| produced n = 4.3t\'^. 



The value a£>*(0) = —1.02 ± 0.05 ~ —1 obtained in our analysis also supports the 
use of the (1 — |xf|) 3 production model which was claimed to describe the available data 



reasonably well |56 



4.5.2 Spectroscopy 

i) Masses of ss states. Our predictions for the masses of pure ss states cannot be directly 
compared to the values for physical states quoted by |28| , because the physical states 



emerge upon mixing with the nonstrange isoscalar states of the corresponding meson 
multiplets. Comparison should be therefore made with the pure ss states calculated from 
both the established physical nn and ss states. Pure ss states were calculated in [p7|] 
with the help of Schwinger's quartic mass formula, for both linear and quadratic masses, 
in essential agreement in both cases and consistent with the quark model motivated 
linear mass relation M(nn) + M(ss) = 2M(sn). Here we quote the masses of the pure 
ss states found in []57| (specifically, from Schwinger's formula for quadratic masses) for 
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vector .1 = 1, 3, tensor .1 = 2 and pseudoscalar .1 = 2 trajectory multiplets, and calculate 
the pure ss mass for the axial-vector J = 1 multiplet as well (all masses are given in 
MeV). We find that these values and the corresponding values in Tables II, IV, VI, VIII are 
in essential agreement: 

• Vector: J = 1 M(ss) = 1014.5 ± 0.4, J = 3 M(ss) = 1862.3 ± 9.0 

• Tensor: J = 2 M(ss) = 1539.4 ± 5.5 

• Pseudoscalar: J = 2 M(ss) = 1869.1 ± 25.3 

• Axial-vector: J = 1 M(ss) = 1487.9 ± 30.0 

Interestingly enough, if one applies the quark model motivated relation mentioned 
above to calculate the mass of the axial- vector sn state using the predicted pure ss mass 
(this work) and the physical a\ (as pure nn) mass (ref. ||28|| ) as inputs, one obtains 
M(sn) = 1359 ± 25 MeV, in good agreement with the corresponding value in Table IV. 

ii) A new cu meson, -D(2637), with mass 2637 ± 2 ± 6 MeV, was recently reported by 
the DELPHI collaboration f58| . Page j59j has demonstrated from heavy quark symmetry 
that the width of .0(2637) claimed by DELPHI is inconsistent with any bound state with 
one charm quark predicted in the Z)(2637) mass region, except possibly D%, Di or 2S 
D, among the three D2 being the best candidate. Interestingly enough, the mass of D2 
predicted in Table VI, 2692 ± 19 MeV, is consistent with the value of the _D(2637) mass 
measured by DELPHI. 

iii) Our predictions for the mass of the B c , B* states are in a good agreement with the 
predictions of different models compiled by Kiselev et al. [^0| which all, including their 



own prediction, lie in the ranges 6281 ± 38 and 6343 ± 26 MeV, respectively. 

The B r mass that we obtained, 6283 ± 79 MeV, is also consistent with the value 



6.40 ± 0.39 ± 0.13 reported by the CDF collaboration pi 



iv) The value of the B 1 meson mass that we obtained, 5692 ± 104 MeV, is in good 



agreement with that recently reported from experiment |62| : 5675 ± 12 ± 4 MeV. 

v) States near the trajectory termination point. As we argued in Sections 2 and 3, any 
nonlinear trajectory of the form (5), including the square-root, can be expected to fit well 
the lowest lying states, and discrepancies are likely to arise only for a few states near the 
trajectory threshold. There are not enough data to confirm or deny this statement, with 
the possible exception of the a^ trajectory. Our prediction for the mass of a 6 , which may be 
the last state on the true 02 trajectory, and is certainly next-to-last on our corresponding 
square- root trajectory (see discusion in Section 4.2), is 10% lower than the observed value. 
This is in contrast to the 04 mass which is lower only by about 1%. The same is true for 
the K\ mass. (Both a 4 and K± lie about 550 MeV below the corresponding thresholds, 
whereas the observed mass of the a 6 is right at our threshold, in contrast to the calculated 
value which is ~ 200 MeV below threshold.) 

Even though the agreement within the errors is acceptable, in our opinion, this raises 
the possibility that near the termination point, the true hadron trajectory grows slower 
than the square-root trajectory. If this possibility is realized, the Kg is likely to lie about 
10% above the value we predict, and it can be expected to be broad, similar to the 
observed a§. 

It should be mentioned that there are other trajectories for which the second state on 
the trajectory is known, e.g. some of the vector or pseudoscalar trajectories. There are 
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two factors which prevent observing the possible general tendency of the square-root form 
to overshoot the true trajectory. First, these states are further away from the threshold 
than those discussed above, in particular, they are about 770 MeV below threshold for 
the vector trajectories, and about 1200 MeV for the pseudoscalars. This means that the 
effect can be expected to be even less than 1%. Second, in both of these cases, we used 
some of the second states on the trajectories as inputs, thus achieving better agreement 
with more of the higher mass states on a trajectory. 

vi) Additivity constraints and mass relations. Our results for masses are in excellent 
agreement with available data. This is not suprising, however, because the additivity 
constraints (26) and (27) are known to lead to high-accuracy higher-power mass relations 
f63fl . Therefore, the numerical spectroscopy results reported here should be viewed solely 
as a confirmation that the particular square-root form of trajectory does not introduce 
any systematic biases. 

4.5.3 A consequence of trajectory thresholds 

If the hadronic Regge trajectories indeed terminate in £, an intriguing possibility that 
might simplify identification of states arises. Even though the maximum £ is very sensitive 
to the specifics of the potential (in particular, its long distance behaviour) and cannot be 
determined at all by fitting the bound states, the thresholds turn out to be insensitive 
to how the potential approaches its asymptotic value and very well determined by the 
bound state spectra. Recall our numerical simulation of the model in Section 2: The 
bound state spectra were well approximated by both extreme cases, the square-root and 
the logarithmic trajectory. These two forms differ drastically in the maximum £ (finite 
vs. infinite), nevertheless, they both predict T within a few percent of the model's true 
value. We therefore conjecture that one can predict "spectroscopy windows" , i.e. range of 
masses for each particular flavor. Even though we cannot obtain stringent error estimates 
on these predictions, it is conceivable that the deviations from the true values are small. 
The most intriguing consequence of this conjecture is, in our opinion, the following: 
According to our calculations, there are no light quarkonia states beyond about 3.2 GeV. 
(The highest threshold for light quarkonia is jT Va = 3.10 ± 0.11 GeV.) Even though the 
charmed states start around 3 GeV, their mixing with glueballs is small. Therefore, if 
our analysis is right, any state above 3 GeV (to, say, 5 GeV) that does not fit into the 
charmed spectra can be expected to be predominantly a glueball or an exotic. 

5 Summary and Conclusions 

Previously [0 we argued that in QCD, the real part of hadronic Regge trajectories 
should acquire curvature and terminate as a consequence of the flux tube breaking due 
to pair creation. In this paper we addressed the issue of what is the specific form of 
the trajectories. We started with the simple potential model of two heavy quarks in 
a potential which is screened at large distances, and we have shown that the parent 
trajectory formed by the bound states of the system can be equally well approximated 
by the square-root and the logarthmic forms. This is not suprising, once the comparison 
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between the potentials which lead to these trajectories, respectively, and the screened 
potential of the heavy quark model is made. Nevertheless, even though the first two are 
indistinguishable, the square-root form is closer to the real trajectory because its real 
part terminates (as does the real part of the heavy quark model), while the logarithmic 
trajectory grows without bound. This growth is specific to the logarithmic form; any other 
trajectory of our nonlinear form (5) with v 7^ has a termination point. This means that 
any v ^ would approximate the true trajectory both quantitatively and qualitatively. 
Different forms can be expected to lead to subtle differences for higher excited states, 
particularly for those near the threshold. Until a distinction can be made, we choose the 
square-root trajectory to study real hadronic spectra. A nice feature of this trajectory 
is that the additivity of inverse slopes reduces to a simple expression which has natural 
units of mass, and that the parameter 7 is given simply in units of inverse mass. It is 
possible that any of these forms, including the square-root, grow faster near the threshold 
than the true trajectory. It is also plausible, although not conclusive, that an indication 
of this tendency has been seen in the 0,2 trajectory. This is also supported by the rate of 
growth of the lattice potential as compared to the potential which leads to the square- 
root trajectory, although the lattice potential should not be taken too seriously at large 
distances. 

With the form of trajectory chosen, the parameters need to be fitted. We have typi- 
cally used masses of few lowest lying states, and/or the intercepts if known and reliable, 
such as the intercept of the p trajectory. We also utilized the additivity requirements. 
Our calculation is in an excellent agreement with various data, both spectroscopic and 
scattering, and it is self-consistent (see the discusion of pseudoscalar b- mesons). 

We conclude that we have provided both strong phenomenological arguments and 
theoretical considerations which indicate that the hadronic Regge trajectories are essen- 
tially nonlinear, and that they can be well approximated by the square-root form (at 
least for the lowest lying states). This observation has a profound effect on our theoret- 
ical understanding of hadronic spectra, namely, it implies that the linear confinement is 
not sufficient and most important factor determining the position of poles, not even for 
states as low as the second or third on a trajectory. Perhaps an even more significant 
consequence of the nonlinearity of trajectories, and, in particular, the existence of the 
trajectory thresholds, is that it may simplify identification of states in experiements. For 
example, once the existence of thresholds is established, it would be easier to identify 
exotic states as those that do not fit in the heavy-quarkonia spectra. 
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A The allowed values of v 

It can be shown that for n < is < n + 1 (for integer is, a(t) = Re oc(t)), the trajectory of 
the form (5) satisfies the following dispertion relation with n + 1 subtractions: 

a {n) (0)t n t n+l r™ , . Ima(i') /A , 

a(t) = a(0) + a'(0)* + ... + — ^- + — jf T dtf ^ 1{ ^ ty (A.l) 

Since for the trajectory (5) Ima(t) = sin(7rz/)(t — T)", the requirement of the positivity 



of the trajectory imaginary part (this requirement follows from unitarity [23| leads to 
2k < v < 2k + 1, where fc is integer. 
For the trajectory parametrization 

a(t) = - 7 (-t)"pog(-t)]0 (A.2) 



(which is the t » T limit of (5), up to a power of logarithm) it was shown in [22| that 
< i/ < 1 and < (3 < 2. 

The range of is can be further restricted. The class of dual models called dual am- 
plitudes with Mandelstahm analycity (DAMA) 0] leaves open a corridor for possible 

asymptotic behavior of a(t), bounded from above by J\t\ and from below by log(|t|) 
J24J1 . (DAMA has the Veneziano limit a(t) ~ t, but the transition to this limit occurs 



discontinuously p^] .) Thus, in DAMA the range of is is squeezed down to < is < 1/2. 
Hence we consider the value of is restricted to lie between and 1/2 p4| . 

B The dynamics of the generalized string model 

By varying the action of the generalized string with massive ends (here the dot stands for 
the derivative with respect to r, Lorentz invariant evolution parameter for the string), 

S gen = ffdrdsL(x,x,x') + J2 L {m) (x t ), (B.l) 

JJ 8=1,2 

one obtains the equations of motion of the generalized string, 

d dL d dL dL 

dr dx ds dx' dx ' 

and the boundary conditions which represent the equations of motion of the massive ends: 

d_dj}^ = dL_ x = ^ 
dr dii dx' 1 

In the gauge r = t discussed above, the equations of motion of the generalized string 
reduce to 

d dL d dL dL . . 

dt dx ds dx' dx ' 
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and the boundary conditions are 

d x\ 

777-1 




d X2 


x\ 


Utf^ 


x 2 



— s = 
dx r 

r)T 

-, . = „. (B.5) 

Let us show that, similarly to the standard case of the string with constant ten- 
sion, there are solutions to the equations of motion of the generalized string (with the 
Lagrangian given in (14)) in the form of a rigid rod connecting the massive ends and 
rotating with frequency u about its center of mass: 

x(t,s) = p(s)(cos(uit), sin (cut), OJ. (B.6) 

Indeed, a = a(\x\) = er(p), since x 2 = p 2 ; therefore dp/dx = x/p = (cos(ui), sm(ut), 0), 
and da/dx = da /dp dp/dx = da /dp (cos(o>£), sin(co>£), 0). Hence 

dL dL da da , r — — / . . . , . \ ._ „. 

to = Tali = - 5? P'V / ^V(-M). -M), 0). (B.7) 

Since also 

d dL auj 2 pd / . . . , . \ ,„ „. 

A*=-7T^?( COSM) ' 8mM) '°)' (B ' 8) 

d <9L ( da . r r— ; r auj 2 pp' 



X^-(--pV 1 -^V+ „ y l, cosM, sinM), , (B.9) 



(the last relation is obtained via da/ds = da /dp p'), it follows that the equations of 
motion (B.4) are satisfied. 

One can show that for the rotation (B.6), the energy of the generalized string is given 
by 

Similarly, the orbital momentum of the generalized string is 

J=J z = jds (x Py - yp.) =jds ™P^ = J fr g _^ . (B.ll) 



Interestingly enough, in his book |64| Perkins also presents the above relations for the 
energy and orbital momentum of the generalized string. He does not however derive 
these relations from the first principles Lagrangian, as in Eq. (14). 

By adding the contribution of the massive ends, one finally has the expressions for the 
total energy and orbital momentum of the generalized string with massive ends: 

B=r-#24,+ r-^4,+ "2_+ "s_, (B.i2) 

JO VI - W 2 p 2 JO VI - U 2 p 2 ^1 - u;2 r 2 ^1 _ uj2 r 2 
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r i dpa(p)up 2 r r2 dpa(p)ujp 2 miur 2 m 2 ujr 2 

o VI - ^ 2 P 2 Jo VI - w 2 p 2 /i _ Cc ,2 r 2 . /i _ w 2 r 2 ' 



Note that the boundary conditions (B.5) define the separations of the massive ends from 
the center of mass through the following nonlinear equations: 



rrijOJ 2 r; 



1 — uj 2 r 2 



a(ri)Jl-cu 2 rl 2 = 1,2. (B.14) 
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